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ABSTRACT 

We report on the first 10 identifications of sources serendipitously detected by the Nuclear Spectroscopic Telescope 
Array ( NuSTAR ) to provide the first sensitive census of the cosmic X-ray background source population at >10 keV. 
We find that these NuSTAR-detected sources are ~ 1 00 times fainter than those previously detected at >10 keV and 
have a broad range in redshift and luminosity (z = 0.020-2.923 and L io— 40keV ^ 4 x 10 41 — 5 x 10 45 erg s '); the 
median redshift and luminosity are z ~ 0.7 and L 10-40 keV ~ 3 x 10 44 erg s , respectively. We characterize these 
sources on the basis of broad-band ~ 0.5 -3 2 keV spectroscopy, optical spectroscopy, and broad-band ultraviolet- 
to-mid-infrared spectral energy distribution analyses. We find that the dominant source population is quasars 
with L 10 — 40 keV > 10 44 erg s 1 , of which ^50% are obscured with Ah > 10 22 cm' 2 . However, none of the 
10 NuSTAR sources are Compton thick (Ah > 10 24 cm' 2 ) and we place a 90% confidence upper limit on the 
fraction of Compton-thick quasars (Lio- 40 keV > 10 44 erg s" 1 ) selected at >10 keV of <33% over the redshift 
range z = 0. 5-1.1. We jointly fitted the rest-frame % 1 0-40 keV data for all of the non-beamed sources with 
10 — 40 kcV > 10 43 erg s' 1 to constrain the average strength of reflection; we find R < 1.4 for T = 1.8, broadly 
consistent with that found for local active galactic nuclei (AGNs) observed at >10 keV. We also constrain the 
host-galaxy masses and find a median stellar mass of ~ 1 0 1 1 M 0 , a factor ~5 times higher than the median stellar 
mass of nearby high-energy selected AGNs, which may be at least partially driven by the order of magnitude 
higher X-ray luminosities of the NuSTAR sources. Within the low source-statistic limitations of our study, our 
results suggest that the overall properties of the NuSTAR sources are broadly similar to those of nearby high-energy 
selected AGNs but scaled up in luminosity and mass. 

Key words: galaxies: active - galaxies: high-redshift - infrared: galaxies - X-rays: diffuse background - 
X-rays: general 
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1. INTRODUCTION 

The cosmic X-ray background (CXB) was first discovered in 
the early 1960s (Giacconi et al. 1962), several years before the 
detection of the cosmic microwave background (CMB ; Penzias 
& Wilson 1965). However, unlike the CMB, which is truly 
diffuse in origin, the CXB is dominated by the emission from 

23 NASA Postdoctoral Program Fellow. 


high-energy distant point sources: active galactic nuclei (AGNs), 
the sites of intense black-hole growth that reside at the centers 
of galaxies (see Brandt & Hasinger 2005; Brandt & Alexander 
2010 for reviews). A key goal of high-energy astrophysics is 
to determine the detailed composition of the CXB in order to 
understand the evolution of AGNs. 

Huge strides in revealing the composition of the CXB 
have been made over the past decade, with sensitive surveys 
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undertaken by the Chandra and XMM-Newton observatories 
(e.g., Alexander et al. 2003a; Hasinger et al. 2007; Brunner 
et al. 2008; Luo et al. 2008; Comastri et al. 2011; Xue et al. 
2011). These surveys are so deep that they have resolved 
**70%-90% of the CXB at energies of %0.5-8 keV (e.g., 
Worsley et al. 2005; Hickox & Markevitch 2006; Lehmer 
et al. 2012; Xue et al. 2012), revealing a plethora of obscured 
and unobscured AGNs out to z ^ 5-6. However, although 
revolutionary, Chandra and XMM-Newton are only sensitive 
to sources detected at ~ 0.5- 10 keV, far from the peak of 
the CXB at ~ 20-30 keV (e.g., Frontera et al. 2007; Ajello 
et al. 2008; Moretti et al. 2009; Ballantyne et al. 2011). Until 
recently, the most powerful observatories with sensitivity at 
~20-30 keV have only resolved ~ 1 %-2% of the CXB at these 
energies (e.g., Krivonos et al. 2007; Ajello et al. 2008, 2012; 
Bottacini et al. 2012) and therefore provide a limited view of 
the dominant source populations (e.g., Sazonov & Revnivtsev 
2004; Markwardt et al. 2005; Bassani et al. 2006; Treister et al. 
2009b; Bird et al. 2010; Tueller et al. 2010; Burlon et al. 2011). 

A great breakthrough in resolving the peak of the CXB is the 
Nuclear Spectroscopic Telescope Array ( NuSTAR ) observatory. 
NuSTAR was successfully launched on 2012 June 13 and is 
the first >10 keV orbiting observatory with focusing optics 
(Harrison et al. 2013). NuSTAR' s focusing optics provide a 

1 order of magnitude improvement in angular resolution and 
a s«2 orders of magnitude improvement in sensitivity over 
previous-generation >10 keV observatories, a revolutionary 
leap forward in performance. One of the primary objectives 
of NuSTAR is to complete a sensitive extragalactic survey 
and identify the source populations that produce the peak of 
the CXB. 

The NuSTAR extragalactic survey comprises three compo- 
nents (see Table 6 of Harrison et al. 2013): a deep small-area 
survey in the Extended Chandra Deep Field-South (E-CDF-S; 
Lehmer et al. 2005) field, a medium wider-area survey in the 
Cosmic Evolution Survey (COSMOS; Scoville et al. 2007) held, 
and a large area (typically shallow) serendipitous survey con- 
ducted in the fields of other NuSTAR targets, including ~ 1 00 
Swift- BAT identified AGNs. In this paper we report on the first 
10 spectroscopically identified sources in the NuSTAR serendip- 
itous survey. In Section 2 we present the NuSTAR observations 
of the serendipitous sources, the multi-wavelength data, and the 
details of our data processing approaches, in Section 3 we de- 
scribe our analysis of the X-ray and multi-wavelength data, in 
Section 4 we present our results, and in Section 5 we outline our 
conclusions. We adopt Hq = 71 km s -1 Mpc” 1 , YIm = 0.27 
and £2 a = 0.73 throughout. 

2. NuSTAR OBSERVATIONS 
AND MULTI- WAVELENGTH DATA 

NuSTAR is the first high-energy (>10 keV) orbiting obser- 
vatory with focusing optics and has a usable energy range of 
3-79 keV (Harrison et al. 2013). NuSTAR consists of two co- 
aligned X-Ray Telescopes (XRTs; focal length of 10.14 m) 
which focus X-ray photons onto two independent shielded focal 
plane modules (FPMs), referred to here as FPMA and FPMB. 
Each FPM consists of four CdZnTe chips and has a ~ 1 2' x 12' 
field of view at 10 keV; the pixel size is 2"46. The focusing optics 
provide NuSTAR with a »1 order of magnitude improvement 
in angular resolution over previous observatories at >10 keV; 
the full width at half-maximum (FWHM) of the point-spread 
function (PSF) is »H8" and the half-power diameter is ~58". 
The absolute astrometric accuracy of NuSTAR is ±5" (90% 


confidence) for bright X-ray sources and the spectral resolution 
is «0.4 keV (FWHM) at 10 keV. 

2.1. The NuSTAR Serendipitous Survey 

The NuSTAR serendipitous survey is the largest-area compo- 
nent of the NuSTAR extragalactic survey program. The serendip- 
itous survey is built up from NuSTAR-detected sources in the 
fields of NuSTAR targets, similar in principle to the serendipitous 
surveys undertaken in the fields of Chandra and XMM-Newton 
sources (e.g., Harrison et al. 2003; Kim et al. 2004; Watson et al. 
2009). A major component of the NuSTAR serendipitous sur- 
vey are ~ 1 5-20 ks observations of ~ 1 00 Swift- BAT identified 
AGNs, which provide both high-quality high-energy constraints 
of local AGNs and «s2-3 deg 2 of areal coverage to search for 
serendipitous sources. However, the serendipitous survey is not 
restricted to these fields and the NuSTAR observations of targets 
not in the E-CDF-S, COSMOS, and Galactic-plane surveys are 
used to search for serendipitous NuSTAR sources; the exposures 
for these targets are also often substantially deeper than the 
NuSTAR observations of the Swift-liXY AGNs (up to on-axis 
exposures of 177.1 ks in the current paper). The expected areal 
coverage of the NuSTAR serendipitous survey in the first two 
years is »G^1 deg 2 . 

Using the NuSTAR data processing and source detection 
approach outlined below, at the time of writing we have 
serendipitously detected »50 sources in the fields of ~70 
NuSTAR targets. Here we present the properties of the first 
10 spectroscopically identified sources; see Table 1. These 
10 sources were selected from NuSTAR observations taken 
up until 2013 January 31. The selection of these sources 
for spectroscopic follow-up observations was based on their 
visibility to ground-based telescopes and they should therefore 
be representative of the overall high-energy source population. 

2.1.1. Data Processing and Source Searching 

The level 1 data products were processed with the 
NuSTAR Data Analysis Software (NuSTARDAS) package (v. 
0.9.0). Event files (level 2 data products) were produced, cal- 
ibrated, and cleaned using standard filtering criteria with the 
nupipeline task and the latest calibration files available in the 
NuSTAR CALDB. The NuSTAR observations of the Geminga 
held comprised 15 separate exposures, which we combined us- 
ing ximage v4.5.1; 24 the other NuSTAR observations reported 
here were individual exposures. 

We produced 3-24 keV, 3-8 keV, and 8-24 keV images us- 
ing dmcopy from the Chandra Interactive Analysis Observa- 
tions (ciao) software (v4.4; Fruscione et al. 2006) for both 
NuSTAR FPMs. 25 We also produced exposure maps in each en- 
ergy band for both FPMs, which take account of the fall in 
the effective area of the mirrors with off-axis angle and are 
normalized to the effective exposure of a source located at the 
aim point. 

We searched for serendipitous sources in all of the six images 
(i.e., the three energy bands for each FPM) using wavdetect 
(Freeman et al. 2002) with an initial false-positive probability 
threshold of 1CU 6 and wavelet scales of 4, 5.66, 8, 1 1.31, and 16 
pixels. To be considered a reliable NuSTAR source we require a 
detection to satisfy at least one of two criteria: (1) to be detected 
in at least one of the three images for both FPMA and FPMB; 


24 See http://heasarc.gsfc.nasa.gov/docs/xanadu/ximage/ximage.html for 
details of ximage. 

25 See http://cxc.harvard.edu/ciao/index.html for details of ciao. 
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or (2) to be detected in at least one of the three images in a 
single FPM but to have a lower-energy X-ray counterpart (e.g., 
detected by Chandra, Swift-XKT, or XMM-Newton). Following 
Section 3.4.1 of Alexander et al. (2003a), we also ran 
wavdetect at a false-positive probability threshold of 1 0 4 
to search the six images (i.e., the three energy bands for each 
FPM) for lower significance counterparts of sources already de- 
tected at a false-positive probability threshold of 1 0 6 in any of 
the three energy bands. 

See Tables 1 and 2 for the details of the X-ray data for the first 
10 spectroscopically identified serendipitous NuSTAR sources. 
All of the NuSTAR sources are detected at >8 keV in at least 
one FPM. 

2.1.2. Source Photometry 

We measured the number of counts for each source at 3-24, 
3-8, and 8-24 keV using either a 30", 45", or 60" radius circular 
aperture centered on the 3-24 keV wavdetect position for 
each FPM; the encircled-energy fractions of these apertures 
are ^0.50, ^0.66, and ~0.77 of the full PSF, respectively, for 
a source at the aim point. The choice of aperture is dictated 
by the brightness of the source and how close it lies to another 
source; see Table 2 for the adopted aperture of each source. 
These measurements provide the gross source counts, which 
we correct for background counts to provide the net source 
counts. To obtain a good sampling of the background counts 
while minimizing the contribution to the background from 
the source counts, we measured the background in source- 
free regions using at least four circular apertures of 45" or 
60" radius at least 90" from the source. The gross source 
counts are corrected for the background counts to give the net 
source counts, rescaling for the different sizes of the source 
and background regions. Errors on the net source counts are 
determined as the square root of the gross source counts. Upper 
limits are calculated when a source is not detected in one 
of the six images or if the net counts are less than the 1 a 
uncertainty; 3a upper limits are calculated as three times the 
square root of the gross source counts. See Table 2 for the source 
photometry. 

2.1.3. Source Fluxes 

The source fluxes are calculated using the net count rates 
(i.e., the net counts divided by the source exposure time) and the 
measured X-ray spectral slope, following a procedure analogous 
to that used in the Chandra Deep Field surveys (e.g., Brandt 
et al. 2001; Alexander et al. 2003a). The X-ray spectral slope 
is determined from the band ratio, which we define here as the 
8-24 keV/3-8 keV count-rate ratio. To convert the band ratio 
into an X-ray spectral slope we usedxsPEC vl2.7.1d (Arnaud 
1996) and the response matrix file (RMF) and ancillary response 
file (ARF) of the detected NuSTAR sources; we produced the 
RMF and ARF following Section 2.1.5. We also used xspec and 
the RMF and ARF to determine the relationship between count 
rate, X-ray spectral slope, and source flux in each of the three 
energy bands: 3-24 keV, 3-8 keV, and 8-24 keV. We calculated 
the source fluxes in the three energy bands using the observed 
count rate and the derived X-ray spectral slope; for the faint 
NuSTAR sources with <100 net counts summed over the two 
FPMs, we set the X-ray spectral slope to T = 1.8, consistent 
with the average X-ray spectral slope of the overall sample (see 
Section 4.3). The source fluxes in each band were then corrected 
to the 100% encircled-energy fraction of the PSF and averaged 
over the two FPMs. 


2.1.4. Source Positions 

To provide the most accurate NuSTAR source positions and 
assist in source matching, we calculated a counts-weighted 
source position. This is determined from the 3-24 ke V net counts 
and the 3-24 keV source position in each FPM. If a source is 
only detected in one FPM at 3-24 keV then the position of the 
source in that FPM is used. 26 

2.1.5. Extraction of the X-Ray Spectral Products 

We extracted the NuSTAR data to be used in the X-ray spec- 
tral fitting analyses. The NuSTAR data were extracted using the 
AwSTAR-developed software nuproducts. nuproducts ex- 
tracts source and background spectra and produces the RMF 
and ARF required to fit the X-ray data; the source and back- 
ground spectra were extracted from each FPM using the same- 
sized apertures and regions as those adopted for the source 
photometry. 

For the serendipitous source in the Geminga field (NuSTAR 
J063358+1742.4) we combined the source and background 
spectrum from each of the 15 observations (see Section 2.1.1) 
to produce a total source and background spectrum. We also 
produced an average ARF file for NuSTAR J063358+1742.4 by 
combining the individual ARF files, weighted by the exposure 
time for each ARF, and we used the RMF produced from the 
first observation when fitting the X-ray data. 

2.2. Lower-energy X-Ray Data 

To extend the X-ray spectral fitting constraints and assist 
in the identification of optical counterparts, we searched for 
<10 keV counterparts for each NuSlA /^-detected source using 
Chandra, Swift-X RT, and XMM-Newton observations. Since the 
NuSTAR serendipitous program targets fields containing well- 
known Galactic and extragalactic targets, they all have lower- 
energy X-ray coverage. However, the only lower-energy X-ray 
data available in the IC 751 field is a short (~2.3 ks) Swift- 
XRT observation in which the serendipitous NuSTAR source 
is detected with only 10 counts by XRT, which is insufficient 
to provide useful <10 keV constraints. For all of the other 
NuSTAR sources there are good-quality <10 keV data and, in 
some cases, there was more than one observation available. 
When selecting suitable lower-energy data we preferentially 
chose contemporaneous observations (i.e., observations taken 
within 1 week of the NuSTAR observations), which was the 
case for three sources in our sample (NuSTAR J032459-0256.1, 
NuSTAR J 121027+3929.1, and NuSTAR J 183443+3237.8). In 
the absence of contemporaneous observations we used existing 
lower-energy data where the 3-8 keV flux agreed to within a 
factor of two with the 3-8 keV flux measured from the NuSTAR 
data; see Section 2.3 for more details. 

2.2.1. Chandra, Swift-XRT, and XMM-Newton Observations 

The archival Chandra observations are analyzed using ciao. 
The data were reprocessed using the chandra_repro pipeline 
to create the new level 2 event file, and the Chandra source 
spectra were extracted from a circular region with a radius 
of ^5"- 10". The background spectra were extracted from 

26 We derive the NuSTAR source name from the counts-weighted NuSTAR 
source position, adjusted to an appropriate level of precision (based on the 
NuSTAR positional accuracy), using the International Astronomical Union 
(IAU) approved naming convention for NuSTAR sources: NuSTAR 
JHHMMSS+DDMM.m, where m is the truncated fraction of an arcminute in 
declination for the arcseconds component. 
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Low-count source and F = 1.8 is used to convert the NuSTAR count rates into fluxes. 

Flux at 3-8 keV measured from the lower-energy X-ray data (either Chandra, Swift-XRT, or XMM -Newton', see Table 1) in units of 10~ 13 erg s -1 cm -2 (see Section 2.2). 
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several source-free regions of ~40" radius, selected at different 
positions around the source to account for local background 
variations. 

The Swift-X RT data are reduced using the HEAsoft (v.6.12) 
pipeline xrtpipeline, which cleans the event files using ap- 
propriate calibration files and extracts the spectra and ancillary 
files for a given source position; 27 the source extraction regions 
had radii of ss20". Since the background in the Swift-XRT ob- 
servations is very low, no background spectra were extracted. 

For the XMM-Newton EPIC data we used the Pipeline Pro- 
cessing System products, which are a collection of standard pro- 
cessed high-quality products generated by the Survey Science 
Center. For our analysis we used the Science Analysis Soft- 
ware (SAS v. 12.0.1), released in 2012 June. 28 After filtering 
the event files for high background intervals, we extracted the 
source spectra from a circular region with a radius of ?»20". The 
corresponding background spectra have been extracted using 
circular source-free regions in the vicinity of the corresponding 
source (ss30"-60" radius regions). Using the SAS tasks rmf gen 
and arfgen we also produced the response matrices for each 
source in each of the three EPIC cameras separately (pn, MOS 1 , 
and MOS2). 

2.3. Counterpart Matching 

To provide reliable source identification we matched the 
NuSTAR sources to the <10 keV and multi- wavelength data; 
see Sections 2.2 and 2.4, and Table 3 for the description of the 
data. We searched for multi-wavelength counterparts within 10" 
of the NuSTAR source positions using on-line source catalogs 
and multi- wavelength images; the latter approach is required 
for faint counterparts or for recent data not yet reported in 
on-line source catalogs. The 10" search radius is motivated 
by the absolute astrometric accuracy of NuSTAR (±5", 90% 
confidence, for bright X-ray sources; Harrison et al. 2013) and 
the low count rates for the majority of our sources. 

A lower-energy X-ray counterpart is found within 10" for 
each of the NuSTAR sources; see Table 2. To provide further 
confidence that the X-ray source is the correct lower-energy 
counterpart to the NuSTAR source, we compared the 3-8 keV 
fluxes of the lower-energy source and the NuSTAR source. We 
selected and extracted the lower-energy X-ray data following 
Section 2.2 and we calculated the 3-8 keV fluxes using a power- 
law model in xspec (the model component is pow in xspec); see 
Table 1 for details of the low-energy X-ray data selected for 
each source. The average source flux was calculated for the 
XMM-Newton data when multiple detectors were used (i.e., pn, 
MOS 1 and MOS 2). In Figure 1 we compare the 3-8 keV 
fluxes from the lower-energy X-ray data to the 3-8 keV flux 
from the NuSTAR data. In all cases the fluxes agree within a 
factor of two, demonstrating that we have selected the correct 
lower-energy X-ray counterpart. 

An optical counterpart is also found within 10" of each NuS- 
TAR source; see Table 3. Given the larger intrinsic uncertainty in 
the NuSTAR source position when compared to the lower-energy 
X-ray source position, we also measured the distance between 
the lower-energy X-ray source position and the optical position. 
An optical counterpart is found within 3" (and the majority lie 
within 1") of the lower-energy X-ray source position for all of 
the sources. See Figure 2 for example multi-wavelength cut-out 
images of NuSTAR J183443+3237.8 in the 3C 382 field. 

27 See http://heasarc.gsfc.nasa.gov/docs/software/lheasoft/ for details of 
HEAsoft. 

28 See http://xmm.esa.int/sas/ for details of the SAS software. 



Figure 1 . Comparison of the measured fluxes in the 3-8 keV band be- 
tween NuSTAR and lower-energy X-ray observations (Chandra, Swift-XKY, or 
XMM-Newton)', see Table 1 for details of the data used for each source. The 
filled circles indicate sources where the lower-energy observations were ob- 
tained within 1 week of the date of the NuSTAR observations and the unfilled 
circles indicate sources where the lower-energy observations were obtained 
more than 1 week after the NuSTAR observations. The solid line indicates 
agreement between the fluxes while the dotted lines indicate a factor of two 
disagreement between the fluxes. 

(A color version of this figure is available in the online journal.) 

2.4. Ultraviolet-Radio Data 

To further characterize the properties of the NuSTAR sources 
we used ultraviolet (UV) to mid-infrared (MIR) data. Table 3 
presents the broad-band UV-MIR photometric properties of 
the NuSTAR sources, primarily obtained from existing, publicly 
available all-sky or large-area surveys, including the Galaxy 
Evolution Explorer (G ALEX', Martin et al. 2005), the Digitized 
Sky Survey (DSS; Minkowski & Abell 1963; Hambly et al. 
2001), the Sloan Digital Sky Survey (SDSS; York et al. 2000), 
the Two Micron All Sky Survey (2MASS; Skrutskie et al. 2006), 
and the Wide-field Infrared Survey Explorer ( WISE ; Wright et al. 
2010). The source photometry is provided in its native format for 
all of the sources. The DSS data, provided for sources outside 
of the SDSS, were obtained from the SuperCOSMOS scans of 
the photographic Schmidt plates (Hambly et al. 2001). As rec- 
ommended by the SuperCOSMOS Sky Survey, all photometric 
uncertainties are set to 0.30 mag for those measurements. Where 
publicly available, we also provide Spitzer photometry from the 
Infrared Array Camera (IRAC; Fazio et al. 2004), obtained from 
the post-basic calibrated data (PBCD) products. To avoid the 
effects of source confusion, photometry was measured in 2" 4 
radius apertures on the O'.' 6 pixel” 1 resampled PBCD mosaics, 
and then corrected to total flux density using aperture correc- 
tions from the IRAC Instrument Handbook (v.2.0.2). 29 Several 
sources were observed during the post-cryogenic Warm Spitzer 
phase, and thus only the two shorter wavelength bandpasses 
from Spitzer-IRAC are available. 


29 See http://irsa.ipac.caltech.edu/data/SPITZER/docs/irac/ for details of 
Spitzer-IRAC. 
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Table 3 

Ultraviolet-to-Mid-infrared Source Properties 
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Figure 2. Example multi-wavelength cut-out images to demonstrate some of the quality of the multi-wavelength data. The object shown here is NuSTAR 
J183443+3237.8 in the 3C 382 field, a faint NuSTAR source. The images are (from top left to bottom right): NuSTAR 3-24 keV FPMA, NuSTAR 3-24 keV 
FPMB, Swift-XKT 0.5—10 keV, DSS B-band, WISE band 1 (Wl; 3.4 // in ), and WISE band 3 (W3; 12 j! m); the NuSTAR images have been smoothed with a 6 pixel 
(14" 8) Gaussian. The circle has a radius of 20" and is centered on the NuSTAR source position. 


In several cases we used photometry from different sources, 
which we list below. For NuSTAR J063358+1742.4 we report 
a /-band non-detection, which is measured from 1.56 ks of 
dithered observations obtained with the Florida Infrared Imag- 
ing Multi-object Spectrograph on the Kitt Peak 2. 1 m telescope. 
The data were obtained on UT 2012 October 17 in photometric 
but 1"6 seeing conditions, and the 3er upper limit was calcu- 
lated in a 2" radius aperture; see Table 3 for more details. For 
NuSTAR J145856-3 135.5 we report the A'-hand magnitude from 
Caccianiga et al. (2008). For NuSTAR J 1 8 1428+3410.8 the op- 
tical photometry comes from imaging reported in Eisenhardt 
et al. (2012), calibrated to the SDSS. The WISE 1 2 p m pho- 
tometry for NuSTAR J 1 8 1428+34 10.8 was measured directly 
from the images as this source does not appear in the WISE All- 
Sky Catalog; we do not provide the shorter wavelength WISE 
photometry for this source as it is superseded by Warm Spitzer 
observations. For NuSTAR J 183443+3237.8 we obtained B-, 
R-, and /-band observations using the Palomar 60 inch tele- 
scope (P60) on UT 2013 March 4 in ~2" seeing; the exposure 
time was 300 s in each band, repeated three times with a 60" 
dither. NuSTAR J 183443+3237.8 was well detected in all three 
bands and the reported photometry in Table 3 was measured in 
4" diameter apertures, which has been corrected for PSF losses. 

We also searched for radio counterparts in the NVSS and 
FIRST VLA surveys (Becker et al. 1995; Condon et al. 1998), 
using a search radius of 30" and 15", respectively. NuSTAR 
J121027+3929.1 was detected in both surveys and has a flux 
of /i. 4 ghz = 18.7 ± 0.7 mJy (in the NVSS survey), which 
corresponds to a rest-frame luminosity density of L | 4 GH/ = 
2.2 x 10 24 W Hz -1 (calculated following Equation (2) of 
Alexander et al. 2003b and assuming a radio spectral slope 
of a = 0.8). With the exception of NuSTAR J01 1042-4604.2, 


all of the other sources had at least NVSS coverage but none 
were detected. The rest-frame luminosity density upper limits 
ranged from L|. 4 ghz < 1.8 x 10 20 W Hz' 1 (for NuSTAR 
J032459-0256.1) to L 1 . 4 GHz < 4.3 x 10 24 WHz' 1 (for NuSTAR 
J 1 15746+6004.9), with the majority of the sources having upper 
limits of L i, 4 ghz < 10 23 -10 24 WHz' 1 . 

2.5. Optical Spectroscopy 

Two of the ten serendipitous sources have existing optical 
spectroscopy: NuSTAR J121027+3929.1 has been previously 
identified as a BL Lac object at z — 0.615 (MS 1207.9+3945; 
e.g., Stocke et al. 1985; Gioia et al. 1990; Morris et al. 
1991) while NuSTAR J 145856-3 135.5 has been previously 
identified as a broad-line AGN (BLAGN) at z = 1 .045 (2XMM 
J145857.0-3 13536; Caccianiga et al. 2008). For the other eight 
serendipitous NuSTAR sources we obtained optical spectroscopy 
at the Palomar, Keck, and Gemini-South telescopes. Table 3 
presents basic information about the observations, including the 
instrument and UT date of the observations and in the Appendix 
we provide specific details for each observation. We processed 
all of the optical spectroscopic data using standard techniques, 
and flux calibrated the spectra using standard stars observed on 
the same nights. 

The optical spectra for the eight newly identified NuSTAR 
sources are shown in Figure 3. Clear multiple broad and/or 
narrow emission lines are detected in six sources, showing 
that the redshift identifications are reliable. However, the op- 
tical counterparts for NuSTAR J 115746+6004.9 and NuSTAR 
J063358+ 1742.4 are comparatively faint and the optical spectra 
are therefore of lower quality when compared to the optical spec- 
tra of the other serendipitous sources. NuSTAR J 1 1 5746+6004.9 
has narrow, spatially extended Lya emission as well as 
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Figure 3. Optical spectra for the eight newly identified serendipitous NuSTAR sources; the optical spectra of the other two sources (NuSTAR J121027+3929.1 and 
NuSTAR J145856-3 135.5) have been previously presented in Morns et al. (1991) and Caccianiga et al. (2008). The prominent emission and absorption lines are 
indicated; see Section 2.5. 


somewhat broadened C hi] emission indicating z = 2.923; spa- 
tially extended Lya emission is often found to be associated with 
powerful AGNs (e.g., Reuland et al. 2003; Geach et al. 2009; 
Yang et al. 2009). The redshift of NuSTAR J063358+1742.4 is 
less certain due to the identification of a single narrow emission 
line, which is more likely to be [Oil] at z = 0.891 than Lya 
due to the rising optical continuum and lack of a strong Lya 
forest decrement (as would be expected had the source been at 
z ~ 4.8); the identification of two absorption features at the 
wavelengths expected for Ca H+K provides additional confi- 
dence for z = 0.891. We consider all of the redshifts to be 
reliable. 

The two NuSTAR sources with existing optical spectroscopy 
(NuSTAR J121027+3929.1; NuSTAR J145856-3 135.5) have 
optical magnitudes consistent with the eight newly identified 
NuSTAR sources and meet our basic requirement for inclusion 
in this paper (i.e., sources identified in NuSTAR observations 
taken up until 2013 January 31); we note that several of the 
other ~40 serendipitously detected NuSTAR sources also have 
existing optical spectroscopy but have been identified in more 
recent NuSTAR observations and so are not included in this 
paper. We therefore believe that the inclusion of these two 
NuSTAR sources does not bias our overall NuSTAR sample. 


3. DATA ANALYSIS 
3.1. X-Ray Spectral Fitting 

To interpret the X-ray data and provide insight into the in- 
trinsic AGN properties of the serendipitous NuSTAR sources 
(e.g., r and Ah) we fitted the X-ray data using physically mo- 
tivated AGN models. We extracted the NuSTAR data follow- 
ing Section 2.1.5 and the lower-energy X-ray data following 
Section 2.2.1. 

For the three sources with >200 counts in each NuS- 
TAR FPM at 3-24 keV (NuSTAR J01 1042-4604.2, NuSTAR 
J1 15912+4232.6, and NuSTAR J121027+3929.1; see Table 2), 
we grouped the NuSTAR data into bins of at least 40 counts bin' 1 
and used x 2 statistics to find the best-fitting model parameter 
solutions. However, the NuSTAR photon statistics were too poor 
to allow for / 2 statistics for the other seven sources, and for the 
X-ray spectral analyses of these sources we fitted the unbinned 
X-ray data using the C-statistic (Cash 1979). The C-statistic is 
calculated on unbinned data and is therefore ideally suited to 
low-count sources (e.g., Nousek & Shue 1989). However, since 
the data need to be fitted without the background subtracted, it is 
essential to accurately characterize the background and use that 
as a fixed model component in the X-ray spectral fitting of the 
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Figure 4. Example X-ray spectra and best-fitting power-law model solutions for NuSTAR J1 15912+4232.6 (left) and NuSTAR J1 21027+3929.1 (right). The NuSTAR 
data are plotted for NuSTAR J1 15912+4232.6 over 4-32 keV and the NuSTAR and Swift-XKT data are plotted for NuSTAR J121027+3929.1 over 0.5-50 keV. The 
black and red crosses are from NuSTAR FPMA and FPMB, respectively, and the green crosses are from Swift-XKT. The X-ray data have been grouped and fitted using 
a power-law model and x 2 statistics; see Section 3.1. The best-fitting models are plotted as solid lines and the bottom panels show the deviations of the data from the 
best-fitting model (A / 2 ). 

(A color version of this figure is available in the online journal.) 


source spectrum. We characterized the background by fitting the 
background regions using a double power-law model (the model 
components are pow*pow in xspec). The photon statistics were 
also often poor for the lower-energy X-ray data (<200 counts) 
and we therefore typically fitted the unbinned <10 keV data 
using the C-statistic with the measured background as a fixed 
component. In Figure 4 we show example NuSTAR spectra for 
two of the brightest NuSTAR sources: NuSTAR J 1 159 12+4232.6 
and NuSTAR J121027+3929.1; for NuSTAR J 121027+3929. 1 
we also show the Swift-XRT data. All fit parameter uncertainties 
are quoted at the 90% confidence level (Avni 1976). 

We initially fitted only the NuSTAR data using a simple power- 
law model (the pow model in xspec) to provide constraints on 
the overall X-ray spectral slope (T) over 4-32 keV. We also 
restricted the NuSTAR data to cover the rest-frame 10-40 keV 
energy range for each source and fitted a power-law model 
to measure both the rest-frame 10-40 keV spectral slope 
(T io— 40 kev) and luminosity (Tio- 40 kev); given the redshift of 
NuSTAR J1 15746+6004.9 (z = 2.923) we fitted to the rest- 
frame 15-60 keV data. See Table 4. 

To provide direct measurements on the presence of absorption 
we jointly fitted an absorbed power-law model (the model 
components arezwABS*POW in xspec) to both the NuSTAR 
and lower-energy X-ray data for each source. 30 For five of 
the sources we fitted the 0.5-32 keV data (for NuSTAR 
J121027+3929.1 we fitted the 0.5-50 keV data, given the good 
photon statistics of this source), jointly fitting the X-ray spectral 
slope and absorbing column density for both of the NuSTAR 
FPMs and the lower-energy X-ray data. However, for NuSTAR 
J1 15912+4232.6 no good-quality low-energy X-ray data exist 
and we therefore only fitted the NuSTAR data, while for the 
remaining three sources (NuSTAR J1 15746+6004.9, NuSTAR 
J 145856-3 135.5, and NuSTAR J 1 8 1428+34 10.8) the photon 
statistics of the NuSTAR data were too poor to provide reliable 
constraints on both T and Ah, and we therefore fitted the 
absorbed power-law model to just the lower-energy X-ray data. 
The best-fitting model parameters are given in Table 4. 

30 We note that AGNs often require more complex models to characterize 
their X-ray emission than that of a simple absorbed power law (e.g.. Winter 
et al. 2009; Vasudevan et al. 2013). Flowever, the data quality of our sources is 
not sufficient to reliably constrain such models on a source by source basis (see 
Section 4.3 for more detailed average constraints). 


3.2. Ultraviolet-Mid-infrcired Spectral Energy 
Distribution Fitting 

To constrain the relative contributions from AGN activity 
and the host galaxy to the UV-MIR data we fitted the broad- 
band UV-M1R spectral energy distributions (SEDs) using the 
0.03-30 jtm empirical low-resolution templates for AGNs and 
galaxies of Assef et al. (2010). Each SED is modeled as the 
best-fit non-negative combination of three galaxy templates and 
an AGN template. The reddening of the AGN template, param- 
eterized by E(B — V ), is a free parameter in the fit. The errors 
on the parameters were calculated using a Monte Carlo method, 
where the photometry is resampled 1000 times according to the 
photometric uncertainties and the SED fits and parameters are 
recalculated; the errors refer to the standard deviation for all of 
the realizations. Since the templates have been empirically de- 
fined using AGNs with similar X-ray luminosities and redshifts 
as the NuSTAR sources, we do not expect there to be significant 
systematic uncertainties in the best-fitting model solutions; the 
efficacy of the SED-fitting approach will be further explored in 
S. M. Chung et al. (in preparation). We refer the reader to Assef 
et al. (2008, 2010, 2013) for further details. 

In Figure 5 we present the UV-MIR SEDs and best-fitting 
solutions and in Table 3 we provide the following best-fitting 
parameters: d (the fractional contribution to the overall emission 
from the AGN component over 0.1-30 jtm; Assef et al. 2013), 
E(B — V) (the dust reddening of the AGN component), L (, /In , 
(the luminosity of the AGN component at rest-frame 6 gm), and 
M : , (the stellar mass of the host galaxy). The stellar mass is cal- 
culated from the absolute magnitude of the stellar component us- 
ing the color-magnitude calibration of Bell et al. (2003). Three 
of the NuSTAR sources have photometric measurements in ^5 
bands (NuSTAR J1 1042+4604.2, NuSTAR J063358+1742.4, 
and NuSTAR J145856-3 135.5) and the derived properties for 
these sources are therefore poorly constrained. 

4. RESULTS 

In analyzing the NuSTAR sources we predominantly fo- 
cus on characterizing their X-ray and UV-MIR proper- 
ties and comparing these properties to those of sources 
detected in previous-generation >10 keV surveys (e.g.. 
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Figure 5. Broad-band UV-MIR SED and best-fitting model solution for the serendipitous NuSTAR sources. The data are fitted with the Assef et al. (2010) AGN 
(magenta dashed curve) and galaxy (elliptical: red dotted curve; spiral: green long-dashed curve; irregular: cyan dash-dotted curve) templates. The best-fitting solution 
is plotted as a black solid curve. The source redshift, best-fitting dust-reddening solution ( E(B — V)), and uncertainties are shown. 

(A color version of this figure is available in the online journal.) 


Swift-BXI', Tueller et al. 2008, 2010; Baumgartner et al. 
2013). 

4.1. Basic Source Properties 

The 8-24 keV fluxes of the NuSTAR sources are up to 
~ 1 00 times fainter than sources previously detected at > 1 0 ke V 
(/8-24keV ~ (0.6-5. 9) x 10~ 13 erg s' 1 cm~ 2 , as compared to 
/ 8 - 20 keV > 0.4 x 10~ n erg s -1 cm~ 2 ; e.g., see Table 2 and the 
RXTE data in Revnivtsev et al. 2004). The NuSTAR sources also 
have fainter optical counterparts and lie at higher redshifts than 
sources previously detected at >10 keV ( R sa 16-22 mag and 
a median redshift of z ^ 0.7, as compared to V ~ 10-16 mag 
and a median redshift of z ^ 0.03; see Beckmann et al. 2009; 
Table 3). 

In Figure 6 we plot the rest-frame 1 0-40 keV luminosity ver- 
sus redshift of the NuSTAR sources and compare them to AGNs 
detected in the Swift-BXI survey (e.g., Burlon et al. 2011). 
With a median luminosity of Z. 10-40 keV ^ 3 x 10 44 erg s -1 , 
the NuSTAR sources are more luminous than the vast majority 


of the Swift-BXI AGNs, where »80% have Z, 10-40 keV < 
10 44 erg s -1 ; the median luminosity of the Swift-BXI AGNs 
is L 10-40 keV ~ 3 x 10 43 erg s -1 . The larger fraction of lumi- 
nous AGNs detected by NuSTAR , in comparison to Swift-BXI, 
is a consequence of the higher sensitivity of NuSTAR and two 
additional factors (1) the strong redshift-dependent evolution 
of luminous AGNs (e.g., Ueda et al. 2003; Barger et al. 2005; 
Hasinger et al. 2005; Aird et al. 2010), and (2) the compara- 
tively small cosmological volume in which NuSTAR is sensitive 
to AGNs with Z, 10-40 kev < 10 44 erg s” 1 (z < 0.2). 

The range of redshifts for the NuSTAR sources is large 
(z = 0.020-2.923). At z = 2.923, NuSTAR J1 15746+6004.9 
is the highest-redshift AGN detected to date at >10 keV that 
does not appear to be strongly beamed (e.g., Beckmann et al. 
2009; Burlon et al. 2011; Malizia et al. 2012). By comparison, 
NuSTAR J032459-0256.1 has a redshift typical of those of 
the Swift- BAT AGNs (z = 0.020) but, with L 10-40 keV ~ 
5 x 10 41 erg s -1 , it is »30 times less luminous than the faintest 
Swift-BAT AGNs; in Section 4.4 we show that this source is 
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Redshift 


Figure 6. Rest-frame 10-40 keV luminosity vs. redshift for the NuSTAR 
sources (circles) compared to the Swift-BAT AGN sample of Burlon et al. 
(2011; squares); filled symbols indicate narrow-line AGNs (NLAGN) and 
open symbols indicate broad-line AGNs (BLAGN). The rest-frame 10 — 40 keV 
luminosity for the NuSTAR sources is calculated directly from the X-ray spectra 
(see Table 4) while the rest-frame 10 — 40 keV luminosity for the Swift - BAT 
AGNs is calculated from the observed-frame 15-55 keV flux, assuming T = 1 .8 
for the /^-correction factor. The dashed lines indicate different flux limits and 
show that the NuSTAR sources are up to ~ 100 times fainter than the Swift-BAT 
AGNs. 

(A color version of this figure is available in the online journal.) 

also unusual since it is hosted in a low-mass dwarf galaxy. 
The high X-ray luminosities for the majority of the NuSTAR 
sources indicate that they are AGNs. However, the origin of 
the modest X-ray luminosity of NuSTAR J032459-0256.1 is 
less clear and it is possible that the X-ray emission is produced 
by a hyper-luminous X-ray source (HLX; e.g., Farrell et al. 
2009; Swartz et al. 201 1) as opposed to a low-luminosity AGN; 
high-spatial resolution observations with Chandra would be 
able to distinguish between an off-nuclear HLX and an AGN or 
nuclear HLX. The median and range in X-ray luminosity and 
redshift of the NuSTAR sources are consistent with expectations 
(Ballantyne et al. 201 1). However, we note that both the redshift 
and X-ray luminosity of NuSTAR J032459-0256.1 are below 
the range typically explored in the models. 

The optical spectral properties of the NuSTAR sources are 
relatively diverse; see Figure 3 and Table 3. Five of the ten 
(^50^22%) serendipitous sources have broad emission lines 
and are classified as BLAGNs, four (^40^%) have narrow 
emission lines and we classify them as narrow-line AGNs 
(NLAGNs), and one is a BL Lac object, with strong power-law 
optical continuum emission and weak emission lines. 31 ’ ’ The 
BL Lac object (NuSTAR J121027+3929.1) is a relatively well 
studied high-frequency peaked BL Lac object (HBL; Padovani 


31 All errors are taken from Tables 1 and 2 of Gehrels (1986) and correspond 
to the ltr level; these were calculated assuming Poisson statistics. 

32 We note that our classification of NLAGNs is fairly loose since we lack the 
emission-line diagnostics around Hey for the majority of our sources to prove 
that they lie in the AGN region of an emission-line diagnostic diagram as 
opposed to the Hu region (e.g., Baldwin et al. 1981; Veilleux & Osterbrock 
1987). 


& Giommi 1995), originally identified at X-ray energies by 
Einstein (MS 1207.9+3945; e.g., Gioia et al. 1990; Morris et al. 
1991 ; Urry et al. 2000; Maselli et al. 2008). Two of the NLAGNs 
have Lkj— 40 keV > 10 44 erg s” 1 and are therefore type 2 quasars, 
representing ~20t 2 1 3% of the NuSTAR sample; by comparison 
six type 2 quasars are identified in the 199 Swift- BAT sample of 
Burlon et al. (201 1), just Rs3^ 2 2 % of the entire sample. However, 
the difference in the fraction of type 2 quasars between NuSTAR 
and Swift-BAT is at least partly related to the increased fraction 
of luminous AGNs in the NuSTAR serendipitous sample; we 
note that, since we lack coverage of the Ha emission line for 
the type 2 quasars, we cannot rule out the presence of broad 
Ha in some of the NuSTAR type 2 quasars. The overall fraction 
of BLAGNs and NLAGNs in the Swift- BAT AGN sample is 
consistent with that found for the NuSTAR serendipitous sample: 
~5(L j % of the Swift- BAT sources are BLAGNs (including 
all Seyfert Is and Seyfert 1.2s) and «5(L 5 5 % are NLAGNs 
(including all Seyfert 1.5s, Seyfert 1.8s, Seyfert 1.9s, and Seyfert 
2s). Therefore, within the limitations of our small sample, the 
biggest differences between the basic properties of the NuSTAR 
sources and the Swift-BAT AGNs appear to be luminosity and 
redshift. 

4.2. X-Ray Spectral Properties: The Presence of Absorption 

The >10 keV sensitivity of NuSTAR allows for the selection 
of AGNs almost irrespective of the presence of absorption, up 
to high absorbing column densities of Ah ^ (1-3) x 10 24 cm -2 . 
However, particularly when using lower-energy X-ray data, we 
can measure the absorbing column densities of the NuSTAR 
sources using the X-ray band ratio (the 8-24 keV to 3-8 keV 
count-rate ratio) and from fitting the X-ray spectra over a broad 
energy range. 

In Figure 7 we show the X-ray band ratio versus redshift 
for the NuSTAR sources and compare them with those expected 
for absorbed power-law emission from an AGN. As can be 
seen, given the high X-ray energies probed by NuSTAR, the 
evidence for absorption can only be clearly identified on the 
basis of the X-ray band ratio for the most heavily obscured 
AGNs (A h > 5 x 10 23 cm' 2 ) at z < 0.5. The X-ray band 
ratios for all of the NuSTAR sources are consistent with Ah < 
5 x 10 23 cm -2 . However, more detailed constraints on the X-ray 
spectral properties and the presence of absorption can be placed 
by directly fitting the X-ray spectra of the NuSTAR sources, 
particularly when including lower-energy data (<3 keV), which 
are more sensitive to column densities of Ah < 10 23 cm -2 . 
We extracted the X-ray spectral products and fitted the X-ray 
data of the NuSTAR sources with an absorbed power-law model 
(zwabs*pow inxsPEC), following Section 3.1; see Footnote 30 
for caveats on the application of an absorbed power-law model 
to characterize AGNs. In Figure 8 we plot the best-fitting X-ray 
spectral slope (r) and absorbing column density (Ah) for the 
NuSTAR sources (see Table 4 for the best-fitting parameters) and 
compare them to the X-ray spectral properties of the Swift-BAT- 
detected AGNs in Burlon et al. (2011). The best-fitting X-ray 
spectral slopes of the NuSTAR sources are broadly consistent 
with those found for well-studied nearby AGNs (T ~ 1.3-2. 3; 
e.g., Nandra & Pounds 1994; Reeves & Turner 2000; Deluit & 
Courvoisier 2003; Piconcelli et al. 2005; Burlon et al. 2011). The 
source with the steepest X-ray spectral slope (T = 2.4 Hq 1 1 5 4 ) 
is NuSTAR J 121027+3929.1, the HBL previously identified at 
<10 keV (e.g., Gioia et al. 1990; Morris et al. 1991). Indeed, 
steep X-ray spectral slopes are typical of HBLs (e.g., Sambruna 
et al. 1996; Fossati et al. 1997). 
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Figure 7. X-ray band ratio vs. redshift for the NuSTAR sources; see Figure 6 
for a description of the symbols. The error bars indicate the la uncertainty on 
the band ratio, which is calculated following the “numerical method” in Section 
1.7.3 of Lyons (1991). The shaded regions show the range of expected band 
ratios for AGNs with Ah < 10 23 cm -2 and Ah « 5 x 10 23 cm -2 , for an 
intrinsic spectral slope of T = 1 .8 ± 0.5; the dotted and solid curves indicate the 
maximum extents in band ratio for Ah < 10 23 cm -2 and Ah = 5 x 10 23 cm -2 , 
respectively. The predicted band ratios were calculated using the RMF and ARF 
for NuSTAR J1 83443+3237.8 (see Section 2.1.5). 

(A color version of this figure is available in the online journal.) 



Figure 8. Best-fitting X-ray spectral parameters (T vs. Ah) for the NuSTAR 
sources and the Swift-BXT AGNs in Burlon et al. (2011). See Figure 6 for a 
description of the symbols; the error bars indicate 90% confidence uncertainties 
for one interesting parameter. The shaded region indicates the range of properties 
found for local AGNs (see Section 4.2). 

(A color version of this figure is available in the online journal.) 


Four of the ten sources (?s40 ^ 9 %) require the presence of 
absorption, with Ah > 10 22 cm -2 , and the other six sources have 
absorbing column density upper limits. The fraction of X-ray 
absorbed AGNs with Ah > 10 22 cm" 2 in the Swift-BAI sample 



Figure 9. Luminosity vs. best-fitting absorbing column density (Ah) for the 
NuSTAR-detected sources and the Swift-BAI AGNs in Burlon et al. (201 1). See 
Figure 6 for a description of the symbols; the error bars indicate 90% confi- 
dence uncertainties for one interesting parameter. The shaded region indicates 
the source properties expected for obscured quasars (L 10-40 keV 10 44 erg s -1 
and Ah ^ 10 22 cm -2 ). 

(A color version of this figure is available in the online journal.) 


of Burlon et al. (2011) is ~53 +4 4 %, indicating no significant 
difference in the fraction of absorbed AGNs between the 
NuSTAR sources and the Swift-BAI AGNs. Eight of the NuSTAR 
sources are quasars with Lio- 40 keV > 10 44 erg s _1 , and four 
(«50t. 4 4%) of the quasars are absorbed with Ah > 10 22 cm -2 ; 
see Figure 9. The fraction of obscured quasars is in broad 
agreement with that found at >10 keV in the local universe 
and from Chandra and XMM-Newton surveys at higher redshift 
(e.g., Ueda et al. 2003; La Franca et al. 2005; Akylas et al. 
2006; Hasinger 2008; Burlon et al. 2011; Malizia et al. 2012); 
however, better source statistics are required to provide sufficient 
constraints to distinguish between different X-ray background 
synthesis models (Gilli et al. 2007). Two of the X-ray absorbed 
quasars are BLAGNs and two are NLAGNs and we discuss the 
origin of the obscuration toward these sources in Section 4.4. 

None of the NuSTAR sources appear to be absorbed by 
Compton-thick material (Ah > 10 24 cm 2 ), despite the 

near obscuration-independent AGN selection over the NuSTAR 
energy range. However, the absorbing column densities of 
Compton-thick AGNs are so high that even the >10 keV emis- 
sion can be significantly absorbed (e.g., AGNs with Ah > 
5 x 10 24 cm 2 can be suppressed by an order of magnitude; 
see Figure 1 1 of Burlon et al. 201 1). Therefore, Compton-thick 
AGNs can be comparatively rare even in high-energy AGN sam- 
ples . 3 ’ Indeed, on the basis of the results obtained for local AGNs 
at >10 keV with the International Gamma-Ray Astrophysics 
Laboratory and Swift-BAI observatories, only ~5%— 10% of the 
detected sources are Compton-thick AGNs (e.g., Tueller et al. 
2008; Beckmann et al. 2009; Burlon et al. 2011; Ajello et al. 

33 Less direct approaches are often required to identify Compton-thick AGNs 
with Ah > 3 x 10 34 cm~ 2 (e.g., optical-mid-infrared spectroscopy, 
photometry, and SED fitting; Risaliti et al. 1999; Alexander et al. 2008; Treister 
et al. 2009a; Goulding et al. 2011; Del Moro et al. 2013; Luo et al. 2013). 
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2012), despite the intrinsic fraction of Compton-thick AGNs 
likely being substantially larger. 34 If distant AGNs have a simi- 
lar range of absorbing column densities as those found locally, 
we would therefore expect ^0.5-1 Compton-thick AGNs in our 
small sample; given the tentative evidence for an increase in the 
fraction of obscured AGNs with redshift (e.g.. La Franca et al. 
2005; Ballantyne et al. 2006; Treister & Urry 2006; Brightman 
& Ueda 2012), we may expect the Compton- thick AGN fraction 
to be even larger in the distant universe. Taking account of the 
low number statistics of our sample, we can therefore place an 
upper limit on the fraction of Compton-thick AGNs of <23% 
in our sample (90% confidence; see Table 1 of Gehrels 1986). 
The 90% upper limit on the fraction of Compton-thick quasars 
over the redshift range of z — 0.5-1. 1 is <33% if we only con- 
sider the seven NuSTAR sources with Z. 10-40 keV > 10 44 erg s _1 . 
These upper limits are marginally too high to distinguish be- 
tween different model predictions for the fraction of Compton- 
thick AGNs detected in NuSTAR surveys for a range of AGN 
luminosity functions and column-density distributions (<22%; 
Ballantyne et al. 2011; Akylas et al. 2012). Better source statis- 
tics are clearly required to accurately measure the fraction of 
distant Compton-thick AGNs. 

4.3. X-Ray Spectral Properties: The Presence of Reflection 

A unique aspect of the NuSTAR data is the insight that it places 
on the > 10 keV emission from distant AGNs and the presence of 
spectral complexity beyond that of simple power-law emission 
(e.g., a reflection component), particularly at z < 1 where the 
rest-frame energy coverage of Chandra and XMM-Newton is 
comparatively modest. By focusing on > 10 keV emission, the 
effect of absorption on the observed emission will be negligible 
(at least up to Ah ^ 5 x 10 23 cm~ 2 ) and the presence of reflection 
can be revealed by the flattening of the intrinsic power-law 
component. 

To investigate the >10 keV emission in our sources we fitted 
the rest-frame 10-40 keV emission using a simple power-law 
model (thepow model inxsPEC), following Section 3.1; see 
Table 4. The spectral constraints for individual sources are poor 
and range from T io- 40 keV ^ 0.4-2. 4, with large uncertainties; 
the mean X-ray spectral slope is rio- 40 keV ^ 1-9. However, 
we can place accurate average spectral constraints by jointly 
fitting the data. When jointly fitting the data we fitted the rest- 
frame 10-40 keV data of the NuSTAR sources with a power-law 
model, jointly fitting the power-law component but leaving the 
normalization for each source to vary independently. In this 
analysis we excluded NuSTAR J 1 2 1027+3929. 1 . the HBL, and 
NuSTAR J032459-0256.1, the low-luminosity system, since we 
wanted to focus on luminous non-beamed AGNs. The best- 
fitting X-ray spectral slope from the joint spectral fitting is 
F io—40keV = 1.88 40 0 26 5 , in good agreement with the intrinsic 
X-ray spectral slope found for nearby AGNs studied at > 10 keV 
(e.g., Deluit & Courvoisier 2003; Dadina 2008; Molina et al. 
2009; Burlon et al. 2011); see Table 5. To first order, the 
comparatively steep average rest-frame 10-40 keV spectral 
slope suggests that there is not a significant reflection component 
in these sources, on average, which would manifest itself as a 
relatively flat X-ray spectral slope at >10 keV (e.g., Nandra & 
Pounds 1994). 

34 Assuming that the intrinsic distribution of absorbing column densities over 
Nn = 10 22 — 10 26 cm' 2 is flat (e.g., Risaliti et al. 1999) and that >10 keV 
surveys are only sensitive to the identification of AGNs with 
Ah < 3 x 10 24 cm' 2 , the intrinsic fraction of Compton-thick AGNs would be 
^20%-40%. 


Table 5 

Joint-fitting Model Parameters 


Model 3 

Sources 13 

Go— 40keV C 

R d 

POW 

PEXRAV 

PEXRAV 

8 

8 

8 

1 qq+0.26 

1 -oo._o.25 
1.8 e 

2.08 + _° 0 25 4 

<1.4 

1.0 e 


Notes. 

3 xspec model used in the joint-fitting process. 
b Number of sources used in the joint-fitting process — the 
low-luminosity system NuSTAR J032459-0256.1 and the 
HBL NuSTAR J121027+3929.1 were not included in the 
joint-fitting process. 

c Best-fitting spectral slope over the rest-frame 10 40 keV 
range. 

d Best-fitting reflection parameter (R: see Footnote 36 for a 
description) over the rest- frame 10—40 keV range. 
e Parameter fixed at given value. 

We can more directly constrain the average strength of the 
reflection component by jointly fitting the rest-frame 10-40 keV 
data using the pexrav model in xspec (Magdziarz & Zdziarski 
1995). 34 Fixing the X-ray spectral slope to r = 1.8 and adopting 
the default parameters for pexrav we constrain the average 
strength of the reflection for the eight NuSTAR sources to be 
R < 1.4. 36 Conversely, if we fix R — 1, the typical value 
found for nearby AGNs selected at >10 keV (e.g., Deluit 
& Courvoisier 2003; Dadina 2008; Beckmann et al. 2009; 
Molina et al. 2009), we constrain the intrinsic X-ray spectral 
slope to be r = 2.08 4 q 2 2 5 4 , also consistent with that of nearby 
AGNs; see Table 5. To first order, our results therefore suggest 
that the strength of reflection in distant luminous AGNs is 
consistent with that found for local AGNs. However, better 
source statistics are required to more accurately constrain the 
strength of a reflection component in distant AGNs and to search 
for changes in the reflection component within subpopulations 
(e.g., dividing the samples in terms of luminosity and absorbing 
column density). 

4.4. Ultraviolet-Mid-infrared Source Properties 

The UV-MIR data of the NuSTAR sources can provide 
insight into the emission from the AGN and host galaxy and 
the presence of dust reddening. Below we first explore the 
MIR colors of the NuSTAR sources and we then analyze their 
UV-MIR SEDs. 

4.4.1. Infrared Color Analysis 

Various work over the past decade has shown that MIR colors 
provide a powerful method to robustly select luminous AGNs 
in a manner that is relatively unbiased by obscuration (e.g.. 
Stern et al. 2005, 2012; Assef et al. 2010, 2013; Donley et al. 
2007, 2012). As such, MIR selection has some similarity to 
hard X-ray selection, and MIR and hard X-ray source selection 
are potentially the two most promising avenues for uncovering 
the full census of AGNs in the universe. Each wavelength has 
various strengths and weaknesses. In particular, various work 
has shown that MIR selection preferentially identifies the most 
luminous AGN with quasar-level luminosities (e.g., Donley 


35 The pexrav model calculates the expected X-ray continuum spectrum due 
to the reflection of power-law emission by neutral material. 

36 The reflection parameter R indicates the solid angle of a neutral slab of 
material illuminated by the primary X-ray source: R w Q./2n. 
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et al. 2007; Eckart et al. 2010), while X-ray selection efficiently 
identifies moderate-high-luminosity AGNs (e.g., Barger et al. 
2003; Szokoly et al. 2004; Xue et al. 2011). On the other hand, 
MIR surveys have now mapped the entire celestial sphere, 
identifying millions of robust AGN candidates. In contrast, 
NuSTAR is unlikely to map more than as 10-20 deg 2 over its 
entire mission lifetime. In order to explore this MIR-X-ray 
complementarity in the new regime offered by NuSTAR, we 
therefore briefly discuss the MIR colors of the 10 serendipitous 
NuSTAR sources. 

Only one of the 1 0 NuSTAR sources (NuSTAR J06335 8+ 1 742.4) 
has four-band Spitzer-IRAC detections, a requirement for the 
Spitzer MIR AGN selection criteria; NuSTAR J063358+1742.4 
is fainter than the WISE flux limits but has IRAC colors that 
place it within the IRAC AGN wedge of Stern et al. (2005). 
Of the other nine NuSTAR sources, eight have at least two- 
band detections by WISE. Stern et al. (2012) and Assef et al. 
(2013) have recently developed WISE AGN selection criteria, 
effectively extending the Spitzer selection criteria across the 
full sky (see also Mateos et al. 2012; Wu et al. 2012). Five of 
the eight NuSTAR sources have WISE colors indicative of an 
AGN according to those criteria. The outliers include the two 
sources with the weakest AGN component (i.e., lowest a val- 
ues; see Section 3.2), NuSTAR J01 1042-4604.2 and NuSTAR 
JJ032459-0256.1. These are the only sources with a < 0.5, 
confirming that MIR selection misses sources where the AGN 
is not bolometrically dominant. 

The final outlier is the HBL NuSTAR J121027+3929.1, 
a BL Lac-type blazar. Massaro et al. (2011) have recently 
published a series of papers discussing the WISE colors of 
blazars. While flat-spectrum radio quasar type blazars have 
colors typical of other AGN populations (e.g., Yan et al. 2013), 
BL Lac-type blazars have unique colors. However, as NuSTAR 
J121027+3929.1 is only detected in the two shorter wavelength 
bandpasses of WISE , it is not possible to compare this source to 
the color criteria developed by Massaro et al. (2011) and Yan 
et al. (2013); note also the caveat in Footnote 3 of Stern & Assef 
(2013). 

4.4.2. Spectral Energy Distribution Analysis 

To quantify the UV-MIR emission of the NuSTAR sources we 
fitted the broad-band SEDs following Section 3.2; see Figure 5 
and Table 3. A significant AGN component ( a > 0.4) is 
required to explain the UV-MIR emission for all of the sources 
except for the low-luminosity system NuSTAR J032459-0256. 1 . 
The rest-frame 6 /im luminosities of the NuSTAR sources 
m ~ (0.9-30) x 10 44 erg s -1 , with the exception of 
NuSTAR J032459-0256.1, which has vL 6/lm « 4x 10 40 erg s" 1 ) 
are in general agreement with that expected for the MIR-X-ray 
(i.e., 6 /rm-2-10 keV) luminosity relationship found for AGNs 
(e.g., Lutz et al. 2004; Fiore et al. 2009); we assumed F = 1 .8 to 
convert between rest-frame 2-10 keV and rest-frame 10-40 keV. 
However, we note that the HBL NuSTAR J 12 1027+3929.1 
and the highest-redshift source NuSTAR J1 15746+6004.9 are 
both X-ray bright compared to the strength of the AGN at 
6 /im, suggesting that the X-ray emission from these sources is 
probably beamed (as would be, at least, expected for an HBL). 

In some cases the presence of dust reddening in the best- 
fitting SED solutions means that the observed contribution of 
the AGN at UV-optical wavelengths is negligible. However, 
we highlight here that, although the strength of the AGN 
continuum at UV-optical wavelengths plotted in Figure 5 is 
inconsistent with the optical spectroscopy in some cases (e.g., 



Figure 10. Dust reddening (E(B — V )) vs. X-ray absorption (Ng) for the 
NuSTAR serendipitous sources from the UV-MIR SED fitting and the X-ray 
spectral fitting, respectively; see Figure 6 for the description of the symbols. 
The solid curve indicates the relationship between dust reddening and X-ray 
absorption expected from the 3 v.Gal — Ml, Gal relationship found in the Galaxy 
(Savage & Mathis 1979; Giiver & Ozel 2009) while the dashed curve indicates 
Av.Gal— 10 x Ml Gal, broadly consistent with that found for AGNs by Maiolino 
et al. (2001). 

(A color version of this figure is available in the online journal.) 

NuSTAR JO 1 1 042-4604.2 and NuSTAR J 1 8 1428+34 1 0.8), they 
are broadly consistent when the range in dust reddening from 
the best- fitting solution is taken into account; see Table 3. As 
expected on the basis of the simplest unified AGN model (e.g., 
Antonucci 1993), the optical emission is heavily extinguished 
in the NLAGNs (E(B — F) ~ 3-6 mag, which corresponds to 
Ay ~ 9-18 mag for Ry = 3.1; e.g., Savage & Mathis 1979), 
with the exception of the low-luminosity system NuSTAR 
J032459-0256.1. There is evidence of dust reddening for two 
of the BLAGNs (NuSTAR J181428+3410.8 has E(B - V) ^ 2 
mag and NuSTAR J 183443+3237.8 has E(B - V) R* 0.6 mag) 
and, as we discuss in the Appendix, the reddening toward 
NuSTAR J1 83443+3237.8 appears to be variable. None of the 
other BLAGNs show evidence for significant obscuration at 
optical wavelengths, as expected for the simplest version of the 
unified AGN model for BLAGNs; however, we note that there is 
a large uncertainty in the dust reddening for NuSTAR J01 1042- 
4604.2, which is due to the limited number of photometric data 
points. 

In Figure 10 we compare the obscuration estimated from the 
UV-MIR SED fitting to that measured from the X-ray spec- 
tral fitting, which provides constraints on the dust-to-gas ra- 
tio in AGNs. The BLAGN with the strongest dust reddening 
(NuSTAR J181428+3410.8) has a measured X-ray absorbing 
column density of /V H ~ 10 22 cm~ 2 , consistent with that 
expected given the A v -V H relationship found in the Galaxy 
(e.g., Giiver & Ozel 2009). The constraints on the X-ray 
absorbing column density for the other four BLAGNs are 
also consistent with that expected given the Ay-Vn relation- 
ship found in the Galaxy; however, in all cases the column- 
density constraints are too weak to rule out the different A v -V H 
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relationship found by Maiolino et al. (2001). By comparison, 
although the three NLAGNs with Lio- 40 keV > 10 43 erg s' 1 
(NuSTAR J063358+1742.4, NuSTAR J1 15912+4232.6, and 
NuSTAR J204021-0056.1) have evidence for significant ob- 
scuration, the inferred X-ray absorbing column densities from 
the dust reddening measurements are lower than those di- 
rectly measured from the X-ray spectral analyses (Nh,a v ^ 
(2 — 4) x 10 22 cm' 2 , as compared to Ah ~ 10 23 cm' 2 measured 
from the X-ray data). However, the dust-to-gas ratios are consis- 
tent with the lower Av-Ah relationship found by Maiolino et al. 
(2001) for AGNs. Neither the HBL NuSTAR J 121027+3929.1 
nor the low-luminosity system NuSTAR J032459-0256. 1 shows 
evidence for significant obscuration in the UV-MIR and X-ray 
bands. 

The best-fitting SED solutions also provide a first-order 
estimate of the host-galaxy stellar masses of the NuSTAR 
sources. The range of stellar masses is large, from ~2 x 
10 9 M q (for the low-luminosity system NuSTAR J032459- 
0256.1) to k» 3 x 10 12 M q (for the highest-redshift source 
NuSTAR J1 15746+6004.9). However, the stellar masses for 
the majority of the NuSTAR sources are relatively tightly 
constrained: the stellar-mass range with these two extreme 
sources removed is (0.7-3. 3) x 10 11 Mq, and the median 
stellar mass is «H0 n Mq. Many of the NuSTAR sources are 
BLAGNs and we caution that reliable stellar-mass constraints 
are challenging for these systems due to the contribution 
of the AGN to the rest-frame optical-near-IR emission (see 
Section 3.2 and Figure 5 for the SED-fitting constraints). 
However, reassuringly, the median stellar mass of the NLAGNs, 
where accurate stellar-mass constraints are less challenging, is 
consistent with that of the BLAGNs when the two extreme 
sources are removed (~ 1 0 1 1 Mq). 

The range and median stellar mass of the NuSTAR sources 
are similar to those of comparably distant AGNs detected at 
<10 keV in Chandra and XMM-Newton surveys (e.g., Babic 
et al. 2007; Alonso-Herrero et al. 2008; Bundy et al. 2008; Xue 
et al. 2010; Lusso et al. 2011). However, by comparison, the 
median stellar mass of the NuSTAR sources is ~5 times higher 
than for z < 0.05 AGNs detected at >10 keV by Swift- BAT 
(«s2 x 10 10 M 0 ; Koss et al. 201 1). To first order this suggests 
that there has been significant evolution in the characteristic 
mass of high-energy emitting AGNs over the redshift range 
Z ~ 0-1 . However, the NuSTAR sources are more luminous than 
the Swift- BAT AGNs and that could bias the results toward more 
massive systems. For example, for a constant average Eddington 
ratio, the order of magnitude higher median X-ray luminosity of 
the NuSTAR sources over the Swift- BAT AGNs (see Section 4.1) 
would lead to an order of magnitude higher black-hole mass and 
thereby a larger stellar mass, assuming no evolution in the black- 
hole-spheroid mass relationship (e.g., Magorrian et al. 1998; 
Marconi & Hunt 2003; Giiltekin et al. 2009). Indeed, Koss et al. 
(201 1) show a weak trend between mean stellar mass and X-ray 
luminosity for the Swift-B.Nl AGNs. Therefore, while our results 
indicate that the most luminous high-energy emitting AGNs at 
Z > 0.1 are hosted by more massive galaxies than high-energy 
emitting AGNs at z < 0.05, a systematic analysis of both local 
and distant AGNs taking account of potential X-ray luminosity 
biases, is required to derive more accurate constraints. 

5. CONCLUSIONS 

We have reported on the first 10 identifications of NuSTAR 
sources serendipitously detected in the extragalactic survey 
program. These NuSTAR sources are ~ 1 00 times fainter than 


AGNs previously detected at >10 keV and have a broad range 
in redshift and luminosity (z = 0.020-2.923 and Lio_ 40 keV ^ 
4 x 1 0 4 1 —5 x 10 45 erg s' 1 ); the median redshift and luminosity 
are z ss 0.7 and Lio- 40 keV ^ 3 x 10 44 erg s' 1 , respectively. 
On the basis of broad-band ~0.5-32 keV spectroscopy, optical 
spectroscopy, and broad-band UV-MIR SED analyses we found 
the following results. 

1. Five (~5(L 34 2 %.j of the ten NuSTAR sources are classified 
as BLAGNs, four («40 432 %) are classified as NLAGNs, 
and one is a BL Lac object. The BLAGNiNLAGN ratio is 
consistent with that found for >10 keV selected AGNs in 
the local universe. See Section 4. 1 . 

2. From fitting the broad-band X-ray spectra we find that the 
dominant source population are quasars with L 10-40 keV > 
10 44 erg s' 1 , of which »50% are obscured with Ah > 
10 22 cm' 2 . However, none of the seven quasars over the 
redshift range z = 0.5-1. 1 are Compton thick and we place 
a 90% confidence upper limit of <33% on the Compton- 
thick quasar fraction. See Section 4.2. 

3. From jointly fitting the rest-frame ~ 1 0-40 keV data for all 
of the non-beamed sources with Lio- 40 keV > 10 43 erg s' 1 
we constrain the high-energy X-ray spectral slope and 
the average strength of a reflection component. We find 
R < 1.4 for r = 1.8 and T = 2.08 4 ° 0 2 2 5 4 for R = 1.0, 
consistent with that found for local AGNs selected at 
>10 keV. See Section 4.3. 

4. From fitting the UV-MIR SEDs we constrain the stellar 
masses of the host galaxies, finding a median stellar mass 
of ~ 1 0 1 1 Mq. The host galaxies of NuSTAR sources are 
~5 times more massive on average than .S'vvz/r-BA'I -detected 
local AGNs at >10 keV. At least part of this implied 
evolution in the characteristic mass of high-energy emitting 
AGNs is likely to be due to X-ray luminosity biases. See 
Section 4.4. 

NuSTAR is providing unique insight into the high-energy 
properties of AGNs, achieving a factor R4 100 times improvement 
in sensitivity over previous observatories at >10 keV. In the 
current study we do not find significant differences in the fraction 
of absorbed AGNs between the NuSTAR sources and nearby 
high-energy emitting AGNs, despite the NuSTAR sources being 
~ 1 0 times more luminous (and ~5 times more massive), on 
average. These results therefore suggest that the central engine 
of distant high-energy emitting AGNs is similar to that of nearby 
AGNs. However, the current study is limited in source statistics 
and provides a first look at the high-energy properties of distant 
AGNs. With the ~ 2 0-40 times improvement in sample size 
afford by the full NuSTAR extragalactic survey (completed in 
the first two years of NuSTAR observations) we will be able 
to make more detailed comparisons and accurately measure 
the high-energy properties of distant AGNs and constrain their 
evolution with redshift. 
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Figure 11. Left: optical spectrum of a Chandra - detected source spectroscopically identified in the Geminga field. The detection of broad Mgn indicates that this 
source is a BLAGN at z = 1.109; right: UV-MIR SED with the best-fitting solution. The data are fitted with the Assef et al. (2010) AGN (magenta dashed curve) 
and elliptical galaxy templates (red dotted curve). The best-fitting solution is plotted as a black solid curve. The source redshift, best-fitting dust-reddening solution 
(E(B — V)), and uncertainties are shown. 

(A color version of this figure is available in the online journal.) 
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APPENDIX 

Here we provide the details of the new optical spectroscopy 
obtained for eight of the serendipitous NuSTAR sources, present 
the optical spectroscopy for an additional Chandra- detected 
source in the Geminga field, and discuss the interesting proper- 
ties of NuSTAR J183443+3237.8. 

A.l. Details of the New Optical Spectrsocopic Observations 

On UT 2012 October 10 we used the Double Spectrograph 
(DBSP) on the Palomar 200 inch telescope to observe NuSTAR 
J183443+3237.8 in the 3C382 field. We integrated for 300 s split 
across two equal exposures in moderate, but non-photometric 
conditions. The observations used the 2"0 wide longslit, the 
6800 A dichroic, the 600/4000 blue grating (e.g., 600 i mm -1 , 
blazed at 4000 A), and the 316/7500 red grating. 


On UT 2012 October 13 we used the DEep Imaging 
Multi-Object Spectrograph (DEIMOS; Faber et al. 2003) at 
the Nasymth focus of the Keck II 10 m telescope to ob- 
serve NuSTAR J20402 1-0056.1 in the AE Aqr field. We ob- 
tained a single 300 s exposure in photometric conditions using 
600/7500 grating. 

On UT 2012 November 9 we used the Low Resolution 
Imaging Spectrometer (LRIS; Oke et al. 1995) at the Cassegrain 
focus of the Keck I telescope to observe NuSTAR J032459- 
0256.1 and NuSTAR J 1 8 1428+34 10.8 in the NGC 1320 and 
W1814+3412 fields, respectively. We observed the sources for 
200 s and 300 s, respectively, in moderate, but non-photometric 
conditions. The observations used the 1"5 wide longslit, the 
5600 A dichroic, the 400/3400 blue grism, and the 400/8500 
red grating. 

On UT 2012 November 20 we again used DBSP at Palomar. 
Conditions were photometric and we used the same instrument 
configuration as employed for NuSTAR J 183443+3237.8 in 
October. We observed NuSTAR J1 15746+6004.9 and NuSTAR 
J1 15912+4232.6 in the SDSS 1157+6003 and IC 751 fields for 
1800 s split into two and three dithered exposures, respectively. 

On UT 2012 December 12 we used the Gemini Multi-Object 
Spectrograph-South (GMOS-S; Hook et al. 2004) at the Gemini- 
South 8 m telescope to observe NuSTAR J01 1042-4604.2 in the 
HLX 1 field. We observed the source for 1200 s, split into two 
exposures dithered by 50 A in central wavelength to fill in the 
chip gap in the focal plane. We used the 1"5 wide longslit and 
600/4610 grating. 

On UT 2013 January 10 we used LRIS at the Keck I telescope 
to observe NuSTAR J063358+ 1742.4 in the Geminga field. We 
observed the source for 1200 s, split into two exposures, using 
the 1"5 wide longslit, the 600/4000 blue grism, the 400/8500 
red grating, and the 5600 A dichroic. The position angle of the 
longslit was set in order to get a second Chandra source in 
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Figure 12. UV-MIR SED and best-fitting solutions for NuSTAR J 183443+3237.8 using (left) our recent (UT 2013 March 4) observations and (right) from the DSS. 
The data are fitted with the Assef et al. (2010) AGN (magenta dashed curve) and galaxy (elliptical: red dotted curve; irregular: cyan dash-dotted curve) templates; the 
best-fitting solution is plotted as a black solid curve. The source redshift, best-fitting dust-reddening solution ( E(B — V)) and uncertainties are shown. 

(A color version of this figure is available in the online journal.) 


the field where there is weak evidence for NuSTAR emission. 
The optical spectrum of this second Geminga serendipitous 
source at a J2 ooo = 06 h 33 m 49*22, <5 J2 ooo = +17°41 , 55'.'l (CXO 
J063349.2+174155) and the UV-MIR SED and the best-fitting 
solution (following Section 3.2) are shown in Figure 11. The 
optical spectrum reveals an AGN at z— 1.109 with somewhat 
broadened Mg n emission, weak [O n] emission, and a strong 
4000 A break with well-detected Ca H+K absorption lines. The 
best-fitting SED solution suggests that the AGN dominates the 
UV-MIR emission. However, since the SED only comprises the 
WISE data, the overall SED is comparatively poorly constrained: 
the best-fitting parameters are a = 0.93 ± 0.03, E(B — V) = 
0.24 ± 0.39, L 6/Mm = (5.34 ± 0.55) x 10 44 erg s' 1 , and 
M* = (9.3 ± 3.2) x 10 10 Mq. There is weak evidence for 
emission from this source in the NuSTAR images. However, 
this source was not formally detected using the source detection 
procedure described in Section 2.1 and we therefore do not 
discuss this source further in this paper. We instead provide 
this information for future researchers of X-ray sources in the 
Geminga field. 

A. 2. Notes on NuSTAR J183443+32378 

NuSTAR J183443+3237.8 is a BLAGN that appears to be 
unabsorbed in the X-ray band (/V H < 1.5 x 10 22 cm -2 ; see 
Table 4). We obtained B- , R-, and /-band observations of this 
field on UT 2013 March 4 using P60; see Section 2.4. The 
optical emission of NuSTAR J 183443+3237.8 has faded since 
the original DSS observations. To explore the origin of this 
fading we fitted the UV-MIR SED of NuSTAR 1 1 83443+3237.8 
following Section 3.2, using both our new data and the older 
DSS data; see Figure 12. On the basis of the original DSS 
observations the best-fitting SED solution indicates E(B—V ) = 
0.00 ± 0.01. However, by the comparison, the best-fitting SED 
solution using the new UV-optical data indicates E(B — V) = 
0.59 ± 0.46, consistent with Ay ~ 1.9 mag for Ry = 3.1 (e.g., 
Savage & Mathis 1979). Assuming the relationship between 
dust reddening and X-ray absorption found in the Galaxy (e.g., 
Giiver & Ozel 2009), the X-ray absorbing column density for 
A v ^ 1.9 mag is Nn « 5 x 10 21 cm~ 2 , a factor «3 below the 


upper limit placed on /V H from the X-ray spectral fitting; see 
Table 4. 
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